On the divergences of inflationary superhorizon perturbations 
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We discuss the infrared divergences that appear to plague cosmological perturbation theory. We 
show that within the stochastic framework they are regulated by eternal inflation so that the theory 
predicts finite fiuctuations. Using the AA'^ formalism to one loop, we demonstrate that the infrared 
modes can be absorbed into additive constants and the coefficients of the diagrammatic expansion 
for the connected parts of two and three-point functions of the curvature perturbation. As a result, 
the use of any infrared cutoff below the scale of eternal infiation is permitted, provided that the 
background fields are appropriately redefined. The natural choice for the infrared cutoff would of 
course be the present horizon; other choices manifest themselves in the running of the correlators. We 
also demonstrate that it is possible to define observables that are renormalization group invariant. 
As an example, we derive a non-perturbative, infrared finite and renormalization point independent 
relation between the two-point correlators of the curvature perturbation for the case of the free 
single field. 



INTRODUCTION 



Primordial perturbations generated during inflation are conveniently characterized by the gauge-invariant curvature 

(— I I perturbation vi^hich can be directly related to the observed CMB anisotropies (]}\. The study of the properties of the 

pi |. curvature perturbation arising in various inflationary models has led to a remarkable confrontation of early universe 

' ' theory with cosmological observations. The increasing observational precision has placed increasing demands on the 

^ , theory and over the past few years results beyond linear theory have been sought, mainly focusing, but not limited 

' to, non-gaussianity. It is well known however that the n-point correlators of the curvature perturbation contain IR 

d ', divergences as a consequence of the approximative scale invariance and gaussianity of inflationary perturbations^. 
Although correlators extending to superhorizon scales as such are not observable, the infrared behavior can manifest 

I [ itself non-trivially both in classical and quantum fleld theory when computing higher order perturbative corrections. 

^ . Moreover, in principle we only have to wait long enough, and the infrared modes will become observable. Therefore, 

' a satisfactory understanding of the theory requires that these IR divergences can be controlled. 

0^ . To demonstrate the problems related to the infrared divergences, it is enough to consider two point correlators and 

g ; single fleld inflation. Using the A., formalism i, the curvature perturbation C(x) can be expressed in a simple form 

H ■ C(x) - iV>(x) + iiV"(<^2(x) - (</>2)) + ... , (1) 
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where TV is the number of e-folds and denotes perturbations of the inflaton fleld on spatially flat sHces. The flrst 
order contribution to the two point correlator of the curvature perturbation behaves as 



>< : (C(xi)C(x2)>(i) = {N'f G(|xi - X2I) ^ -(7V')^7'^ln(^^^-^) , (2) 



where G(|xi — X2I) = (0(xi)0(x2)) is the correlation function of the scalar fleld, V gives the amplitude of the inflaton 
perturbations 4>, and we have considered the case of scale invariant fluctuations. The scale L is an a priori arbitrary 
cutoff needed to regulate the infrared divergences of the two point correlator. In reality, the theory may contain a 
physical scale beyond which scale invariance is broken and L may be related to that scale. For instance, one could 
argue that in our local patch inflation has lasted only a flnite time so that beyond some infrared scale, the assumption 
of scale invariance necessarily breaks down. In any case, in the actual experiments only differences of correlators are 
measured and they remain flnite even in the limit L — > 00 , 

(C(xi)C(x2)>(i) - (C(xi)C(x3))(i) ~ -{N'r-P^ln . (3) 

Xl X3 

The value of the IR cutoff L is therefore irrelevant when considering observable quantities calculated using flrst order 
perturbation theory. 



^ More accurately, the expressions depend on an IR cutoff L. 
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The situation becomes more complicated beyond the Hnear order since the divergent part of the correlators will no 
longer be simply an additive constant. For example, the second order contribution to the two point correlator behaves 
as 

which in the limit L ^ 00 contains divergences depending also on the separation of the two points Xi and X2 . 

There are two distinct issues associated with the appearance of these IR divergences. The first is whether 
G(|xi — X2I) is actually IR divergent. The divergence is the outcome of the assumption that the spectrum is scale 
invariant, an assumption that will probably break down after some scale L. One might expect that the complete 
theory of inflation predicts a well deflned value for G(|xi — X2I), maybe large but free of divergences. The other 
issue is the relation of this (large) value of Gdxi — X2I) to the measurements of an observer accessing only a small 
patch of the whole universe. In other words: how do the predictions of the theory deflned on the largest possible 
scales translate into predictions for measurements on much smaller scales? These problems have been discussed in 
the literature 0, B SSBl! but there seems to be no agreement on the correct approach. 

In this paper we address both issues. First, we show that the variance of the fleld (0^) is flnite and well deflned 
for generic slow-roll inflation. We prove this using the framework of stochastic inflation which accurately describes 
the fluctuations of quantum flelds on superhorizon scales. Furthermore, we note that since IR fluctuations are only 
observationally relevant in regions which have exited the eternally inflating regime and have eventually thermalized, 
the scale of eternal inflation V{(p-Ei) ^ e(0Ei)Afp provides an ultimate upper limit for observable fluctuations. Only 
below the length scale L^i, corresponding to i^ei, does the notion of a well deflned evolving background fleld apply 
and the spacetime can be approximated by a global FRW metric. We then proceed to consider how predictions at this 
scale can be related to predictions for observations restricted in patches of size M < L < Lei by applying a systematic 
renormalization prescription. In particular, we focus on the renormalization of the classical one loop expressions for 
the 2-point and 3-point functions of the curvature perturbation. We flnd that the renormalization amounts to shifting 
the background flelds in the connected parts of the correlators. Furthermore, we note that apart from some contrived 
models, the exact position of the patch of size M in the larger universe is irrelevant in single fleld inflation, to the 
extent that the background evolution corresponding to the observable patch M has been deflned. Therefore, averaging 
over all possible embeddings of M in the whole universe has no meaning in this context. 

The paper is organized as follows: In sections II and III we use the formalism of stochastic inflation to demonstrate 
that for potentials supporting slow-roll inflation, the one point functions of the scalar fleld and the curvature per- 
turbation are flnite as long as the energies are sub-planckian. Therefore there are no real IR divergences associated 
with inflationary fluctuations. Then, in section IV we relate predictions at scales L to scales M < L using the AA^ 
formalism, and argue that a renormalization group equation can be used to relate predictions obtained using two 
different cutoffs/renormalization scales. In section V we discuss a concrete example, the free fleld, and demonstrate 
that it is possible to construct an observable that is both infrared flnite and renormalization point independent. In 
section VI we close with a discussion and a comparison of our flndings to recent work on the subject. 



II. ETERNAL INFLATION AND ONE-POINT CORRELATION FUNCTIONS OF THE INFLATON 

In this Section, we will show how effects of eternal inflation lead to the regularization of the infrared divergences 
present in the correlation functions of the inflaton. Although we will focus on the chaotic inflation scenario [1, our 
analysis can be easily extended to the case of new inflation [1, inflation driven by multiple flelds [lH or to the 
case of the curvaton scenario 

Let us suppose that primordial inflation is driven by a single scalar fleld (j> with the potential V{4>) satisfying the 
slow roll initial conditions e = iM|, iV4,/Vf < 1, 77 = M'^V^^/V < 1. During one Hubble time At ^ the 
value of the inflaton fleld changes by 

On the other hand, at the same time scale fluctuations of the inflaton Scj) are generated with a characteristic wavelength 
I ~ ^ H~^. These fluctuations have a randomly distributed amplitude among different causally disconnected 
regions (Hubble patches), and the width of this distribution is given by 
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As a result, observers living in diflerent Hubble patches see different expectation values of the inflaton, because for 
them long wavelength (superhorizon) fluctuations (50 are physically indistinguishable from the inflaton zero mode. 

In some Hubble patches the sign of the infrared fluctuation Scf) is positive, in others it is negative and therefore may 
overcome the effect of the classical force |[5]) acting on the inflaton. When < \54)\, stochastic fluctuation 5(j) may 
lead to the effective growth of the of expectation value of the inflaton (/) in a given Hubble volume. This regime is 
denoted as eternal inflation because there will be always Hubble patches where the expectation value of the inflaton 
continues to grow |[13t]. 

One can see that the stochastic force becomes more important than the classical one when |(Mp/T^)F0| < ^JV/Mp 



or 

V{M > e(<^Ei)M|, (7) 

where e is the slow roll parameter. For example, in the chaotic inflation model the condition ^ is satisfled when 
Since the upper bound on possible values of field corresponds to the Planckian energy density, 
we conclude that the regime of eternal infiation is realized at A~^/^Mp < < X~^^'^Mp, a very wide interval of 
possible values of the inflaton, since A ^ 10^^'^. 

Larger values of fleld correspond to larger scales of perturbations in the chaotic inflationary scenario, and one may 
conclude that the regime of eternal inflation can influence the infrared behavior of the correlation functions of both 
(p and the curvature perturbation (. In fact, as Woodard has shown leading infrared logarithmic divergences in 
the correlation functions of (p and ( can be reproduced in the formalism of stochastic inflation which is thus well 
suited for treating the regime of eternal inflation. 

In this formaHsm one decomposes the Heisenberg operator for the field (j) into the infrared and ultraviolet parts 
according to the prescription 

${N, x) = $(7V, x) + I <fk 6{k - 5aH) (akMN)e''''^ + al0^(7V)e*'^^) + 5<jy, (8) 

where N = log a is the number of infiationary e-folds, (5 ^ 1 is a small parameter^, 6{x) is the Heaviside step 
function and 5(j) is a small contribution suppressed by additional powers of the slow roll parameters. We are especially 
interested in the dynamics of the infrared part $ of the inflaton fleld, which is responsible for the infrared divergences. 

Substituting the decomposition ([8|) into the equation of motion for the Heisenberg operator (j), one can show that 
the infrared part $ satisfles the Langevin equation 

where / is the composite operator containing contributions of the ultraviolet modes and having the correlation property 

{f{N)f{N'))^^6{N~N'), (10) 

as can be checked by a direct calculation of this correlator in the Bunch-Davies vacuum state. The remarkable 
property of the equation ^ is that all its terms commute with each other and therefore can be considered as classical 
stochastic quantities. The Langevin equation ^ describes the process of the inflaton random walk due to the effect 
of generated superhorizon fluctuations. 

At the next step, by the Stratonovich prescription one proceeds from the Langevin equation ^ to the Fokker-Planck 
equation 

dP _ 1 92 Q f^.dV 



dN 37rAf2 902 (^^) + 8^ 90 (^^ 90 ^j' ^^^^ 

which governs the evolution of the probability P(0, N) to measure a given expectation value of the inflation in a given 
Hubble patch. Its solution can be expressed in the form 

+00 



P(0, N) = «(0)e-''(^) V c„V«(0)e-^"^, (12) 



n=0 



^ This parameter is intrinsically related to the decoherence scale Id, where perturbations which are leaving the horizon finally become 
classical (i.e., to the scale, where decaying mode can be neglected as compared to the growing one). The scale decoherence 1^ corresponds 
to few e-folds after perturbations leave the horizon. Note that 5 cannot be made arbitrarily small as it is related to the slow roll parameters 
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where 



while V'n and En are the eigenfunctions and eigenvalues of the Schrodinger equation given by 



2V'n+2(-« +(1- )')^« =^nV^«^- (14) 



Due to the supersymmetric form of the potential in the Schrodinger equation (|T4|) all eigenvalues satisfy the condition 
En > 0, and it is trivial to see that the ground state corresponds to a zero eigenvalue. To be precise, the ground state 
corresponds to i?o = if and only if the corresponding eigenmode is normalizable. However, this is always the case 
for chaotic inflationary models with V{4>) bounded from below and growing as ^ 0. 

The Fokker-Planck equation (fTTj) should be supplemented with initial conditions (in particular, they define constants 
Cn with n > 1 in the solution (fT2|) ). It is physically reasonable to choose them near the hypersurface of Planckian 
energy density l^((/)max) — Mp corresponding to the cosmological singularity.^ If the universe started its evolution 
from the single Hubble patch of the nearly Planckian size (with, say, Hf ^ ^Mp ), one takes 

mA^O='^('/'-^0max), (15) 

where Ni = log Ui and Oi is the initial value of the scale factor in that Hubble patch. Observe that one cannot choose 
initial conditions of the form Pi = 5(0 — 0max) due to breakdown of the quasiclassical approximation used for the 
derivation of the Langevin equation ^ . 

From the solution lfT2|) we see that in the deep infrared limit N — log a cx3 or, more precisely, at log a '3> E^^ only 
the term with Eq ^ survives in the solution (fT2|) . i.e., the distribution function P{(f>, N) reaches its time- independent 
asymptotics 

PLv(0)^m"'cxp(^-^) (16) 

which is called the Linde-Vilenkin wavefunction of the Universe [l5l.fl^. Therefore, one-point correlation functions of 
the infiaton field are given in the limit — ^ oo by the expression 



d#"FLv(</'), (17) 



where 0min is given by the condition e(0inin) ~ 1, the normalization is defined as N = d(f)P{(f>), and n is arbitrary 

(limited only by the condition that the corresponding integrals converge). 

Therefore, we have explicitly shown that the naively diverging infrared parts of all one-point correlation functions 
of the infiaton field are rendered finite once one takes the effects of eternal infiation into account. 

One may observe that, although being finite, the correlators ((/)")ir still remain very large. As an example, one 
can take chaotic infiationary model with the potential V{(j)) — \m?(fp' . Using Eqs. lfT6|) . (fl7|) and taking into account 
the the upper bound for the value of the infiaton field is 0max — "/^M'^/m, we find that ((/)^)ir ~ 0.96Mp/m^ and 
(</'^)iR ~ 0.15Mp/m* . The level of non-gaussianity generated by eternal infiation ((/)'*)ir — 3((/i^)jj^ « — 2.61Af|,/m^ 
is extremely large. However, this amount of non-gaussianity will never be accessible for an observer within a given 
Hubble patch since the horizon of this Hubble patch will never cross the hypersurface of eternal infiation [l^ . 

For completeness, let us discuss what happens with the Fokker-Planck probability distribution if eternal infiation 
comes to an end before the stationary asymptotics l|16p is reached. 

At < 0EI one has to neglect the first, stochastic, term in the right-hand side of the Fokker-Planck equation (flTI) . 
and the latter acquires the form 

dP Ml d fl dV \ 

-^{t7^p)- (18) 



ON Stt d4> \V d(t) 



^ Contrary to the general lore, the regime of eternal inflation does not solve the problem of cosmological singularity, and the latter remains 
unavoidable even in eternally inflating space-time fl7| . By choosing such initial conditions for the Fokker-Planck equation, we neglect 
Hubble patches where the singularity is achieved but the overall energy density in the patch is dominated by matter flelds other than 
the infiaton field 4>- 
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Its solution is given by 



Pi^,N)^^flN+ f d^^v(^] ), (19) 



where the function / should be determined from the evolution of the probability distribution P in the regime of 
eternal inflation. The physical meaning of the solution lfT9|) is that the probability is transported unchanged along 
the hypersurface 



Note that we count the number of e-folds starting from the beginning of inflation instead of from its end, i.e., we 
define N = log {a/ai), where Qi is the initial value of the scale factor. 

Finally, we note that the multipoint correlation functions of the form {4>{xi)4>{x2) ■ ■ ■ 0(x„)), where xi, X2, ■ • •, Xn 
are separated by distances larger than the coarse graining scale (i.e., the horizon scale) depend on the Starobinsky's 
decoherence parameter S and cannot be calculated without considering carefully what happens at near-horizon scale. 



III. ONE-POINT CORRELATION FUNCTIONS OF THE CURVATURE PERTURBATION 

Let us now show that one-point correlation function of the curvature perturbation (C") also remains finite if one 
takes properly into account the effects of eternal infiation. Our goal will be to calculate the average duration of 
the inflationary stage (A'^totai) as well as its higher order correlation functions (-/Vtotai) arbitrary n. Correlation 
functions of the curvature perturbation at the end of inflation (i.e., at = A'^totai) are trivially related to the latter: 

(C"(^total)) = ((A^total - (A^total))"). (21) 

To calculate these one-point correlation functions, we will use the method developed by Starobinsky in [§| and expand 
it for the case of chaotic infationary scenario. 

It is basically guaranteed that inflation will come to the end in a given Hubble patch when the inflaton expectation 
value drops below the boundary deflned by the Eq. ([7]) , and the evolution of the inflaton fleld becomes determin- 
istic. However, there are still stochastic fluctuations of the inflaton fleld which may lead to a sudden change of in a 
given Hubble patch from a value (f) > 0ei to the value 0min < 4> ^ 0ei- Thus, to flnd the total expected number of 
e-folds for the given Hubble patch one has to use the formalism of stochastic inflation. 

The probability distribution for the end of inflation or, in other words, the total number of e-folds iVtotai can be 
determined from the probability distribution P(0, N) from conservation of probability by 



W^(iVtotal) = P(0, N) 



M' 



2 



— — lim 



P(0,^total). (22) 



Instead of dealing with this distribution function directly, we will calculate the correlation functions 

r+oo 

Q„(0)= / d^^"P(0,7V), (23) 

where Ni = log ai and ai is the value of the scale factor at the beginning of inflation. The moments Q„ can 
be naturally related to the correlation functions of the total number of e-folds (^totai) ^^"^ curvature perturbation 
correlation functions (C") , as we will see later. 

The zero moment Qoiff') can be determined from the equation 

1 ^\ivQ,H^4,WZiyQ.}) = -pi*.^.)- m 



After taking into account that (5o(0 = '/'max) = 0, i.e., that the probability flow through the hypersurface of Planckian 
energy density is absent, one flnds 

^^('/')Qo(0) = 3A/|.e-2-W) r°''^d</)'e2''W') r P((/., M), (25) 
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where v{4>) is again given by Eq. (flSl) . From this expression one can immediately see that the probabiUty distribution 
function l(22|) for the total duration of inflationary stage is properly normalized: 



1 dV 






JNi 



dNw{N) = 1. 



In turn, the higher moments Qn are governed by the recursive set of equations 

1 92 



3A/2 902^' ' Stt 90 ^'^ 



nQn-i, 



(26) 



(27) 



relating higher moments Qn with the lower ones Qn-i- The general solution of Eq. (|27|) can be written in the form 
similar to (l25l): 



^^('/>)Qn(0) - 3Af 



(28) 



We immediately conclude that all the higher moments Qn{4>) are finite, since Qo{(f>) and all the integrals in this 
expression are well behaved. 

The correlation functions of the total number of e-folds are related with moments Qn [(f)) according to the prescription 



/l,f-2 



j_dV_ 
V~d6 



Qn = n 



d4>Qn-i{4>), 



(29) 



where we have used Eq. In particular, the average duration of the inflationary stage is given by 



total/ 



d^Qo{^)^iMl 



(30) 



P(</.,iVO 



(31) 



Finally, the curvature perturbation correlation functions are related to the moments Qn{4>) according to the pre- 
scription 



(r(iVtotal)) = n / dcj)Qn^m 



dcl^Qoi^) 



(32) 



As we see, all one-point correlation functions of A^totai and ^ (computed at the end of inflationary stage) are flnite 
for all physically interesting chaotic inflationary models. The flnal answer does depend on the initial conditions for 
eternal inflation due to the dependence of the zero moment Qo on the initial distribution function P{<j),Ni). For 
example, one has 



total/ 



(33) 



juf2n 

{C{N,o,,,))^Dn^ (34) 

for the m2</)2 chaotic inflation model, where C„, Dn = Cn — Cq ^ C(l) are numerical constants depending on the 
initial conditions. By virtue of the definition Di — 0. 

These large non-gaussianities however never become observable for any given observer within a given Hubble patch, 
and we now turn to a discussion of the observable quantities. 



■* In the deterministic inflationary regime the first term on the right hand side can be neglected with respect to the second one (see the 
end of the previous Section). 
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IV. THE TWO-POINT & THREE-POINT FUNCTIONS OF THE CURVATURE PERTURBATION AT 

ONE LOOP 



In the previous section we demonstrated that the one-point functions of the scalar field and the curvature pertur- 
bation are finite when the effects of eternal infiation are taken into account. Thus, there are no real infinities in the 
theory. Furthermore, the scales alluded to in the previous section are much larger than any currently observable scale. 
Then the question arises: How do the predictions at scales L translate into predictions at scales M < L7 We address 
the issue of relating two different scales in the present section using the prescriptions of the AN formalism [2]. 

According to the AA^ formalism, the perturbative expansion for the correlator of the curvature perturbation to one 
loop is given by [l^ 

(C(xi)C(x2)) = {^NaN^ + NAN^''{cl)^)fl) Gi2 + ^NabN^'^GI^ , (35) 

where Gy = G{rij) = G{\'x.i — Xj|) = ((/)(xi)(/)(xj)). The parameters L and I in the symbol ((?!>^)('2-) are the infrared 
and ultraviolet cutoffs, respectively, appHed to all quantities in equation (|35l) . The 3-point function to one loop is 
given by 



(C(X1)C(X2)C(X3)) = 



NabN^N'' + f -Ni^r'N^'N^ + Nir.N^'^Nr '-^'^^'^ 



(G12G13 + G21G23 4- G31G32) 



v^sjv jv -r i^-i^ABC^'^ ^^cjKVI^L) 

+ -^NabcN^^N'~^(^ (G12 + G13) G12G13 + (G21 + G23) G21G23 + (G31 + G32) G31G3 

+ NabN^'^N^ G12G13G23 , (36) 

where both of the above formulae assume that 

(<^-^(x.)0^(x,))='5^^G.,. (37) 

Thus, G(r) is the correlation function of a generic light field in quasi de-Sitter. The exact scalar field correlator G(r) 
is supposed to be known up to a superlarge scale L and is presumably well defined - see previous section for its 
coincident limit. According to the AA^ prescription, all the derivatives of N in l(35|) . (|36l) are to be evaluated at the 
value of the background field corresponding to some time after the shortest scale of interest has left the horizon. Note 
that the expansion l(35|) requires the existence of a well defined background evolution. Then, the picture of stochastic 
infiation imposes that the largest field value for which the expansion (|35|) is applicable is given by the hypersurface 
of eternal infiation ([7]) with (fiEi- For larger field values, (p > 4>ei, the notion of the background field breaks down and 
so does the expansion ((35|) . Spacetime above this scale can no longer be described by a global FRW metric and the 
correlators G(r) are no longer translationally invariant. However, it is important to note that an observer will never 
interact with the surface of eternal infiation since such an interaction presupposes a universe which has thermaHzed 
as only in the thermalized patches will the horizon grow, permitting such measurements. 

Consider now an observer located in some random position. Since the state of the field is translationally invariant, 
the stochastic properties of the field will be described by the correlation function G(r) for any such observer. However, 
the field fiuctuations will not be experimentally accessible beyond some length scale AI defining the causal patch of 
the observer. Fluctuations with wavelengths longer than M should appear as part of the background for this patch. 
Since these fiuctuations only contribute a constant to G(r) for r < M, the correlator appropriate for describing 
measurements within the observer's patch is 

G(r)^G(r)-(</)2)j^5). (38) 

In other words, an observer limited to measurements in a volume r < M will be able to measure G(r) up to a constant. 
One has 

r/M l/l 
r/L 1/M 

where 'P^{k) is the power-spectrum obtained from the fourier transform of G(r). The first term in G l(39|) gives a 
small finite correction to the correlation function which goes to zero as r becomes much smaller than M . The minimal 
length scale / corresponds to the shortest cosmological scale of interest; it is known to be at least 60 e-folds smaller 
than the size of the causal patch M (or larger depending on the temperature at reheating) . 
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Let us now express Eq. l(35|) in terms of G{r), the correlator appropriate for the patch of the size M. Keeping only 
one loop terms, we have 



(C(xi)C(x2)) = NaN- 



-{NaN'')b 



(L) 



1 



1 Gu + ^-NabN^'^GU + NaN^{^')[^,^ . (40) 



The last term is a constant, which cancels out when comparing differences of two correlators. The coefficient of G{r) 
now includes the term 



1 



^ +:,iNAN^) 

<t>L ^ 



A\B 
B 



(41) 



where (j)L denotes the background field appropriate for the scale L and (/)(x)|^^'' are fiuctuations with wavelengths 
larger than M. The 3-point function gives a similar result. Expressed in terms of G(r) it takes the form 



(C(X1)C(X2)C(X3)) = 



(x)[i? 



abcN^'n'' + N2bN"'-Nc 

BC at a 



tBC-, 



+ -NabcN^^'n'' 



Gi2 + Gi3 ) G12G13 



G12G 



13 



NabN N(j G12G13G23 



+ 2NabN'^N 



Aj.rB/>2\{M) 



G 



12 



G 



13 



G 



23 



const 



(42) 



where we again have retained only terms up to one loop. We thus see that the result of L ^ M is the replacement of 
the tree level coefficients by their averages. 



An{M) ^{An[4>L + 



(43) 



where Ai = N^Na , A2 — NabN^N^ etc, as well as the appearance of a disconnected part which for the two point 
correlator is simply a constant. Otherwise, the functional form of the expansion as given by the rules of [l^ remains 
the same with all integrals now truncated at the new IR cutoff M . The same should hold for all coefficients of all 
correlators in the AA^ expansion and therefore we write 



(C(xi)C(x2))m - {{NaN^) + (iVAA^r>('/'')[5^^)) G{r) + \{NABN^'')G\r) + 



(44) 



and 



(C(xi)C(x2)C(x3))a/ = 



{NabN^'N") + ( -{N^BcN^'N'^) + {N^bN'^'^Nc) ] {4>') 



2\(0 



--{NabcN^''n'^ 



G12G13 



G12 + Gi3 ) G12G13 + ...)+ {NabN^'^NS) G12G13G 



23 



+ 2{NABN^N^){(f')\^L^ (Gi2 + Gi3 + G23) + 



(45) 



Thus, the coefficients appropriate for the cutoff M can be calculated by replacing (pL ^ (f^L + '/'(x)(^'^' and averaging 

over all field fiuctuations (/>(x)|^^^ with wavelengths longer than M. This appears to be the generic form for the AiV 
expansion when used in patches smaller than Lei- 

Equation l(43l) could also be interpreted in the following way: for an observer in the patch of the size M making 
observations at r < M, the quantity 0l + 't'i'^)^^^) ^® ^ slowly varying function, indistinguishable from a background 
field. Of course, this background differs according to the location of the patch in the larger volume of size M. Thus, 
the average in equation (|43l) can be seen as an average over the position of the patch of the size M in the larger 
volume or, equivalently, an average over all possible "backgrounds" seen by observers residing in M @|. Of course, 
this quantity is not (necessarily) related to observable quantities of interest to any particular observer in a specific 
patch, in agreement with the results of 0,0]. 

The above considerations might seem to suggest that the only appropriate infrared cutoff for observers like us is the 
current horizon L — I/Hq 0|. Defining the theory at L = ^/Hq, seems to avoid the need of using averages in the AiV 
expansion. However, framed in this way, the choice of the IR cutoff seems dictated by convenience rather than the 
theory itself and is hence ad hoc, designed to hide the long wavelength fiuctuations. However, fiuctuations on scales 
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larger than our current horizon presumably exist and we "simply" have to wait a few billion years to access them. In 
the future we will be forced to calculate using the cutoff appropriate for that time and thus obtain a different result 
for our local curvature perturbation. Therefore, using a larger box should not introduce any theoretical uncertainty 
for predictions of the local curvature perturbation field. 

Equation (|40l) allows us to deal with a change in the IR cutoff. We can define the coefficients in the expansion l(35|) 
to be evaluated at a background field appropriate for our current patch and impose 1 /Hq as an IR cutoff. If now we 
are granted access to longer wavelengths, say by waiting long enough and repeating our observations, it is expected 
that we will infer a different background evolution since we will now be averaging over a different volume and more 
long wavelength fiuctuations will be visible. Thus, we must accordingly revise the infiationary history that we match 
to the post-infiationary universe. This will of course lead to a variation in the coefficients of the AN expansion for 
the connected part of the correlators, given by 



^^^'h-^-^^B (46) 



This running shows how the connected parts of the correlation functions of the curvature perturbation change as the 
cutoff scale M is increased. It is akin to the renormaHzation group fiow in quantum field theory, where the coupHng 
constants run with the renormalization scale, and we will accordingly refer to M as the renormaHzation scale for 
cosmological perturbations. We note, however, that at least for the 2-point function (|40|) the renormalization simply 
amounts to a constant change in the numerical value of the connected part of the correlator (see (|35l) ). The difference 
of the 2-point functions, which is a measurable quantity, thus remains invariant under the change of the cutoff scale 

(C(xi)C(x2)) - (C(xi)C(x3)) = (C(xi)C(x2))m - (C(xi)C(x3))m . (47) 

Equation l(46|) admits the following interpretation, at least in the case of single field infiation. The coefficients An 
are functions of the derivatives of the number of e-folds N{(j)) w.r.t. 0: An{N', N" , . . . , iV^")). We assume that from 
(|46l) one can derive a consistent set of equations for the running of the derivatives of the form 

|^=^(7V',7V",...,iV(«)). (48) 

We see that as M is changed, the values of the derivatives A^*^"' ((/)), and thus the function iV((/)) itself, are changed 
according to l|48p . Furthermore, equation (|48l) takes a particularly simple form for infiationary models with a monomial 
potential V{(j)) — Kip". Indeed, for such models 

n Mp n Mp 



dN' 1 7V"2 

= -^7^^(VM). (50) 



and eqs (|48|) simply become 



d\nM 2{N 

We illustrate this point with a simple example in the next section and refer to future work for further investigation. 

V. A CONCRETE EXAMPLE: FREE FIELD 

In this section, we consider two point correlators of the curvature perturbation in single field infiation with quadratic 
potential V = ^rri^cfp', although the following discussion would proceed practically unchanged for any monomial 
potential. We have 

All - An 
" = -Ml* ■ " = -Ml ■ 

and all the higher order derivatives vanish. The two point correlator (|35l) is thus written to all orders in loops as 

(C(xi)C(x2)) = {N' f + kN'rCj, . (52) 
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According to l(48|) . the value of N' will change as the cutofl is increased. We parameterize the variation as 

47r - 

N' = -j^^* + f3{M) , P{l/Ho) = , (53) 

where we have assumed that the original cutoff has been set at our present horizon Mq = I/Hq, such that expressions 
([ST]) hold for this cutoff, and 4>* denotes the background field appropriate for this cutoff. We thus have 



df3 1 / 47r y V{l/M) 



d\nM 2\MlJ ^^^+I3{M) 



(54) 



from which we obtain 



W^T^-^^fl- Wl-:5?ln^oM) (55) 



where we have neglected the scale dependence of V . Thus, with the new cutoff M, N' reads 



At: - V 47r - 

iV' = -^0jl-^lni/oA/.-^0M. (56) 



In this case, changing M simply amounts to evaluating expressions ([511) at a new field value, shifted by 



^^4>^l~^lnHoM. (57) 

Note that the form of the theory remains unchanged, ie still described by a quadratic potential. This is a general 
feature of all monomial potentials. 

So far, we have shown that one is free to choose the IR cutoff for infiationary perturbations, provided appropriate 
shifts in the parameters of the theory are performed, see l|57p . It is also interesting to note that for a free field 
it is simple to construct manifestly infrared finite expressions containing the one-point correlators of the curvature 
perturbation. The expression for the two point correlator (|52l) can be inverted to yield 



Gi2 = r (^-1 + ^1 + (87r2)-i(Mp/0)4(C(xi)C(x2)) j , (58) 

where we have neglected the other unphysical root that would yield perturbations of the same order as the background 
field G{r) (fr. The correlator G{r) in (|58|) of course depends logarithmically on the IR cutoff but the difference 

i/i i/i 

( sin(fcr) sin(fcr')\ _ f 'P^{k) A . ,2 



l/L 1/L 

contains only negative powers of L and is thus finite even in the limit L — > 00. Thus we can construct a manifestly 
infrared finite expression 

G12 - Gi3 = ^2 (^^l + (87r2)-l(Mp/^)4(C(xi)C(x2)) - Y^l + (8^2)-l(^Jp/^)4(^(x^)^(x3)) ^ . (60) 

By approximating the field correlator as 

G(r) « -T'^ln (^) , (61) 

one further obtains 

^Jl + (8^2)-i(Mp/0)4(C(xi)C(x2)) - ^1 + (87r2)-i(Mp/^)4(C(xi)C(x3)) = ^In (^^) , (62) 

where the dependence on the cutoff has completely disappeared. This can be seen as a non-perturbative generalization 
of the first order result ([3]), which is rederived by expanding (|62l) to first order in (Mp/(/))^(CC). Since equation i(62|) 
does not depend on the cutoff L at all, it remains unchanged in the renormalization prescription described above and 
is thus renormalization group invariant to the precision of the approximation l(6T1) . Given a measurement of the two 
point correlator of the curvature perturbation (C(xi)C(x2)) for one separation of points ri2, equation l(62|) can thus 
immediately be solved to find an unambiguous prediction for the correlator at any other separation ri3. 
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VI. DISCUSSION 

The appearance of IR divergences in cosmological perturbation theory has been addressed by many authors with 
conclusions that do not always coincide 0, 0, SHH]- For example, Lyth [3| has suggested that the appropriate IR 
cutoff for inflationary calculations is a box slightly larger than the observable universe. He argued that the use of a 
much larger box leads to problems since one would further need to average over the position of our observable patch 
in this larger box and that this would introduce theoretical uncertainties. However, framed in this way, the choice of 
the IR cutoff seems dictated by convenience rather than by the theory itself and is hence ad hoc. On the other hand, 
the authors of 01 have pointed out that IR effects are indeed irrelevant for the observations but conceded that using 
a box larger than the horizon would require one to average over the position of our patch leading to a prediction for 
an averaged, and therefore inherently statistical, power spectrum. In contrast, we do not find any such theoretical 
uncertainty, nor is there any need to average over the possible embeddings of our observable universe into the box 
defined by the hypersurface of eternal infiation, since the field correlator is manifestly translation invariant. This due 
to the fact that the box itself is defined by the end of infiation, i.e. each point has experienced the same fixed number 
of e-folds. 

Using the framework of stochastic inflation, we flrst showed that for slow roll inflation the fluctuations are in fact 
finite to all scales up to Planckian energy density^. However, fluctuations on such scales are unobservable, and only 
regions which have thermalized are interesting observationally. We then related the predictions for the curvature 
perturbation on the largest possible thermaHzed scales to observations performed in much smaller regions by an 
appropriate redeflnition of the coefficients in the AA^ expansion. As we have discussed, this procedure can be carried 
out to the case of 3-point functions, and we conjecture that this is so for any n-point correlator. 

Our flndings show that it is indeed permissible to use our horizon as a cutoff but that there is no uncertainty related 
to using a larger box. In perturbation theory the difference between the two choices is an additive constant plus a 
redeflnition of the background at which the coefficients of the AA^ expansion l(35|) should be evaluated. In a sense, 
the size of the box M represents simply the renormalization point, the redeflnitions a renormalization prescription. 
The largest possible box is deflned by the requirement that the fleld has dropped below the eternal inflation threshold 
(see Eq. ([Tj)) and has therefore thermaHzed. Since G{r) is translationally invariant, any observer conflned in a patch 
of size M is able to probe it for r < M and obtain the same answer regardless of where the patch is located. 

Thus, in general the correlators of the curvature perturbation were seen to depend on the renormalization scale AI. 
Since the (perturbative) renormalization of the two-point correlator amounts to a constant change in the numerical 
value of the connected part of the correlator, the difference of the two-point functions actually remains invariant under 
the change of the cutoff scale. As discussed in section V, in the free single field case it is straightforward to find also a 
non-perturbative generalization of the first order result which is independent of the infrared cutoff L. As is apparent 
from Eq. (|62|) . given a measurement of the correlation of the curvature perturbation between the points xi and X2, 
one can find an unambiguous and renormalization point independent prediction for any (C(xi)C(x3)). It would be 
very interesting to generalize this result to multifield cases as well as to non-trivial potentials. 
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